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Abstract. 

Earlier observations of NGC 1808 in various wavebands (X-ray, optical, near-infrared, radio) provided evidence for 
the existence of either a starburst or a Seyfert 2 nucleus. We here present the results of multiwavelength XMM- 
Newton and Chandra observations, which directly prove the co-existence of thermal diffuse plasma and non-nuclear 
unresolved point-like sources associated with the starburst activity, along with a Low Luminosity Active Galactic 
Nucleus (LLAGN) or an Ultra Luminous X-ray source (ULX). The broad bandwidth of XMM-Newton allows us 
to show that the unresolved nuclear source in NGC 1808 dominates the hard X-ray spectrum, while the emission 
in the soft regime, below 1 keV, is dominated by a thermal component associated to an extended starburst. 
Both EPIC and RGS data provide reliable detections of a number of emission lines from heavy elements, with 
abundances ranging from roughly 0.7 to 2.2 Zq for different elements. However, no 6.4 keV Fe Ka fluorescence 
line emission was detected. The analysis of the nuclear region of NGC 1808 allows us to detect and disentangle the 
contribution of an unresolved nuclear X-ray source and the starburst region, but the exact nature of the nucleus 
remains unknown. The observed luminosity of NGC 1808 is L2-10 keV = (1.61±0.06) x 10 40 erg s _1 . A comparison 
of our OM 212 nm image with a CTIO 4-m telescope Ha frame shows a good general correspondence between the 
emission from massive stars and warm ionized gas, with minor deviations near the ends of the bar in NGC 1808. 
An aditional, very soft thermal spectral component with kT~0.1 keV has been discovered in the XMM-Newton 
spectral analysis, which most likely originates from the halo of NGC 1808. 

Key words. Galaxies: active - Galaxies: nuclei - Galaxies: general - X-rays -ISM: general - Galaxies: evolution - 
Galaxies: halos - Galaxies: starburst - individual: NGC1808 



1. Introduction 

NGC 1808 (a(J2000)=05:07:42.34; <5(J2000)=-37:30:47.0 
with an uncertainty of l'/25, obtained from 2MASS 1 ), 
classified as an SABb Seyfert 2 galaxy, is located at a 
distance of 10.9 Mpc (H = 75 km s" 1 Mpc" 1 , 1"= 53 
pc). Images in different wavebands suggest a high star 
formation (SF) intensity in the central region with a di- 
ameter of ~ 20" (i.e., w 1 kpc). There are several optical 
hot spots associated with H Il-regions (Sersic & Pastoriza 
1965). The luminous and compact knots detected in radio 
and infra-red (IR) images (Saikia et al. 1990, Kotilainen 
et al. 1996) that do not coincide with the optical hot spots 
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are probably supernova remnants (SNRs) or complexes of 
unresolved SNRs. Dust filaments were found perpendic- 
ular to the central disk plane and can be explained as 
outflowing material driven by supernovae (SNe; Heckman 
et al. 1990). A recent interaction of NGC 1808 with its 
companion NGC 1792 could explain both the intense star- 
formation activity and its peculiar morphology (Dahlem 
et al. 1990, Koribalski et al. 1993). 

The nature of the nucleus is still unclear. It is classified 
as Seyfert 2 based on the optical emission lines detected 
by Veron-Cetty & Veron (1985). However, there is also 
evidence from optical observations of a hidden starburst 
(SB), based on polarization measurements by Scarrott et 
al. (1993). The nucleus is claimed to have a strong non- 
stellar component, because only 10% of the IR radiation 
observed by ISO could be interpreted as emission related 
to star formation (Siebenmorgen et al. 2001). On the other 
hand, soft X-ray (0.1-2.4 keV) data obtained with ROSAT 
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(Dahlcm et al. 1994 and Junkes et al. 1995) favor an inter- 
pretation in terms of SF-related emission, but do not dis- 
card other hypotheses. Awaki & Koyama (1993) interpret 
their Ginga X-ray data (1.5-37 keV) as an obscured ac- 
tive galactic nucleus (AGN). However, Awaki et al. (1996), 
using ASCA (2-10 keV) observations, point out that the 
hard X-ray spectrum could also be the result of starburst 
activity, while the long-term variability from the Ginga 
and ASCA observations again suggests a Seyfert nature 
of the nucleus. 

With the advent of the XMM-Newton and Chandra 
satellites, the coexistence of starburst and AGN activity 
has been detected in X-rays in the nuclei of several nearby 
galaxies. NGC 4303 presents evidence of the composite 
SB / AGN nature of the nuclear region: Chandra and UV- 
HST observations show that a superstellar cluster coexists 
with a hard X-ray source, possibly an AGN, in the 3 cen- 
tral pc of the galaxy (Jimenez-Bailon et al. 2003, Colina et 
al. 2002). At larger scales, Persic et al. (2004a) claim that 
starburst and AGN coexist in NGC 4666 by detecting evi- 
dence of diffuse thermal emission, high-luminousity X-ray 
sources, and the presence of an AGN in the X-ray band. 
The detailed analysis of our XMM-Newton data shown 
here, combined with other observations, sheds new light 
on the nature of the nuclear activity in NGC 1808 and the 
associated outflow of gas. 

The observations and data reduction are described in 
Sect. our results are presented in Sect. |31 followed by 
our interpretation of the inner region and circumnuclear 
emission in Sect. l4.il and of the more extended X-ray emis- 
sion outside the nuclear area in Sect. 14.21 The conclusions 
of this work are given in Sect. GD 

2. Observations and data reduction 

2.1. XMM-Newton data 

Information on the XMM-Newton instrumentation is pro- 
vided by Jansen et al. (2001; XMM-Newton mission), 
Striider et al. (2001; EPIC-pn), Turner et al. (2001; EPIC- 
MOS), den Herder et al. (2001; RGS), and Mason et 
al. (2001; OM). Our XMM-Newton observation (Obs- 
Id 0110980801) was performed on April 6, 2002. The 
scheduled exposure time was 40 ks, while effective expo- 
sure times per instrument are listed below. The EPIC-pn 
exposure was in extended full frame mode and with the 
thin filter. Both RGS cameras were in spectroscopic mode 
and the two MOS cameras in full frame, with the thin fil- 
ter. OM was used in full frame low resolution mode with 
the U, UVW1 (219 nm) and UVW2 (212 nm) broad band 
filters, plus the two (optical and UV) grisms. 

The data were processed with the Science Analysis 
Subsystem, SAS, v.5.4.1 (Gabriel et al. 2004) and the most 
up-to-date calibration available in February 2003. The 
standard tasks to process data, ep/mchain and rgsproc, 
were run with default parameters. For the EPIC in- 
struments, checks against pile-up were performed with 
the epatplot task. No sign of pile-up was found for any 



of EPIC-pn, MOS1, or MOS2 data. Additionally, back- 
ground flaring time intervals, for which the signal-to-noise 
does not improve by including these events, were filtered 
following the method described in Piconcclli et al. (2004b). 
Only single and double events were considered in the anal- 
ysis. For RGS the selection was for periods with less than 
0.12 c/s in the background region of the CCD number 9. 
The effective exposure times after each filter was applied 
were 31.4 ks, 39.0 ks, 38.6 ks, and 37.4 ks for EPIC-pn, the 
two MOS cameras, and RGS, respectively. The OM broad 
band data were processed with the SAS omichain task, 
which performed all necessary corrections. Unfortunately, 
about half of the UVW1 filter image was lost due to 
ground-station problems, and the UV grism data was also 
lost due to telemetry drops. 

2.2. Chandra data 

The Chandra ACIS (AXAF CCD Imaging Camera) 
archival image, sequence number 700451, was extracted 
from the Chandra archive. The ACIS observation took 
place on December 19, 2002 with an exposure time of 
43 ks. We used the science level 2 files, generated by the 
pipeline standard processing on December 19, 2002, as 
stored in the archive. They were processed with Ciao 3.0.2. 
According to the information available on the Chandra 
pages, further corrections should be applied (i.e. reduc- 
tion of tap- ringing distortions). However, for the purpose 
of the present work, this was found to be unnecessary. No 
additional corrections were therefore applied. 

2.3. Ground-based Ha imagery 

The Ha data used here were obtained with the CTIO 4-m 
telescope on November 18, 1993. A 2048 x 2048 chip with a 
0'.'431 pixel scale was used, leading to a total field-of-view 
of 147 x 147. 

Three exposures of 20 min. each were taken with a red- 
shifted Ha filter (A = 657.8 nm; AA = 0.64 nm), with 
a total integration time of 1 hr. The continuum subtrac- 
tion was performed using an i?-band image with a total 
integration time of 8 min. The data are extremely sensi- 
tive to low surface brightness emission, but as they are 
uncalibrated, will be used only for qualitative studies. 

3. Results 

3.1. X-ray imaging 

XMM-Newton EPIC X-ray images of NGC 1808 were pro- 
duced by grouping the photons from the cleaned event 
lists into the sub-bands defined by the XMM-Newton 
Survey Science Centre (SSC; Watson et al. 2001) and 
using the merge task in SAS v.5.4.1, which takes ex- 
posure time corrections into account. These images are 
displayed in Fig. ^ The maximum emission is located 
at a EPIC (J2000)=05:07:42.40; c5 EPIC (J2000)=-37:30:46.2 
with an uncertainty of l'.'5-3". The position of the 2MASS 
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Fig. 1. EPIC images of NGC 1808 in the following energy bands (from left to right and from top to bottom): 0.2-0.5 
keV, 0.5-2.0 keV, 2.0-4.5 keV, 4.5-7.5 keV, 7.5-12 keV, and the total bandpass 0.2-12 keV. All images are displayed 
on a logarithmic scale with a threshold of 3d of the mean background value. 



nucleus is indicated by a cross. An adaptive filter was 
used to smooth the images in areas of low signal-to-noise 
(S/N) ratio, while not smearing out the emission from 
point sources. 

NGC 1808 shows extended emission in the soft X-ray 
images up to 4.5 keV (Fig. [I}. Fig.|21shows the comparison 
of the radial profiles of NGC 1808 and MCG-06-30-15 in 
two energy bands: 0.2-4.5 keV and 4.5-10 keV. The QSO 
MCG-06-30-15 can be considered as a point-like source 
in all the XMM-Newton energy bands (Kirsch et al. 2005). 
The comparison shows that the emission of NGC 1808 in 
the lowest energy band is extended while its radial profile 
in the 4.5-10 keV band does not differ from the QSO pro- 
file, taking the errors into account. Therefore, NGC 1808 
can be considered as a point-like source above 4.5 keV. 
This points out that the origin of the X-ray emission is not 
due to the compact central source alone. The elongated X- 
ray emission of NGC 1808 follows the same orientation as 
the Ha emission (see Sect. l3.4|l . The extended morphology 
is detected at energies below 4.5 keV, with a maximum in 
the 0.5-2.0 keV band (see Fig.^J. At these low energies, 
thermal emission processes tracing high-mass SF play an 
important role. Thus, our data suggest that the stellar ac- 
tivity in the nuclear region is significantly contributing to 
the observed soft X-ray emission. With increasing energy, 
the central emission maximum becomes more prominent 
and the non-nuclear sources fade away. 

In order to resolve more details, we have compared 
our 0.2-12 keV EPIC image with the Chandra ACIS data. 



Both were smoothed with an adaptive filter and are dis- 
played in Fig. [3J The Chandra ACIS image shows how the 
nucleus is surrounded by a population of discrete sources. 
The image (left panel of Fig. |3J) shows the detailed struc- 
ture in the central part of NGC 1808: at least four point- 
like sources are detected in the image, apart from the dif- 
fuse emission. Two of them are located at ~4" (212 pc) 
NW and SE from the maximum of the X-ray emission 
in the Chandra image and have also been detected in a 
Chandra HRC observation. Zezas et al. (2001) reported a 
luminosity for each of them of around 10 39 erg s _1 . 



Although the nucleus is observed as a point-like 
source in the HRC image (see Zezas el al. 2001), the ACIS 
image reveals the presence of a double-peak structure 
(see Fig. 0). The hard (E>2 keV) and soft (E<1.5 keV) 
X-ray images of the two peaks presented in Fig. 0] reveal 
that the eastern source, SI, is significantly softer than 
the western one, S2. The locations of the two peaks 
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and o£ handra (J2000)=05:07:42.19; 
^ handra (J2000)= -37:30:45.8 with an uncertainty of 
0'.'6. The nuclear location from 2MASS is indicated by an 
arrow in the Chandra image, closer to the eastern source 
found by Chandra, SI, but compatible within the errors 
with both of them. The remaining two point-like sources 
are situated on the ring of radio emission as imaged 
by Saikia et al. (1990), close to the detected SNR and 
H II-regions. 
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Fig. 2. Comparison between the radial profiles of NGC 1808 and the point-like source MCG-06-30-15 in the 0.2- 
4.5 keV and 4.5-10 keV bands. Y-axis shows the normalised brightness (normalised counts/s/px 2 ). The plots show 
that the emission of NGC 1808 is extended in the 0.2-4.5 keV band and point-like in the 4.5-10 keV band. 




Fig. 3. NGC 1808 0.2-12 keV X-ray images taken with the Chandra ACIS (left panel) and the XMM-Newton EPIC 
(right panel) with the ACIS image overlaid as contours (contour levels are 2,5,10,15,20, and 50 counts). The nuclear 
location from 2MASS is indicated by an arrow in the Chandra image. Both images are displayed on a logarithmic scale 
with a threshold of 3er of the mean background value. 



3.2. X-ray spectral analysis 

We performed an extensive X-ray spectral analysis of 
NGC 1808 using the EPIC and RGS instruments on board 
XMM-Newton, as well as, with Chandra's ACIS camera. 

The XMM-Newton EPIC data were used to perform 
a spectral analysis of NGC 1808 in the 0.35-9 keV en- 
ergy band with moderate resolution (80 eV at IkeV). In 
order to search for differences in the X-ray spectrum of 
the nucleus and its surroundings, we defined three cir- 
cular regions and one annulus. The circles have radii of 
16" (850pc), 35" (1.9kpc), and 2' (6.4kpc), while the 
annulus has external and internal radii of 20" and 50" 
(1.0-2.7 kpc), respectively (excluding one extra-nuclear 
source). Background regions were selected to be on the 
same CCD as the source but far enough to prevent source 
contamination. Figure |5] shows the spectra of each region. 
We simultaneously analysed the EPIC-pn and RGS spec- 
tra, and finally we also studied the Chandra spectra of 



the two nuclear point-like sources SI and S2 (see Fig. 
both located within the 16" inner region extracted in the 
XMM-Newton analysis. 

In order to apply the modified \ 2 minimization tech- 
nique in the spectral analysis, all EPIC spectra were 
grouped such that each spectral bin contains at least 50 
counts. The Chandra spectra were grouped with at least 
20 counts in each bin. The spectra were analysed using 
XSPEC v.11.3.0 (Arnaud 1996). 

The spectra of the 35" and 2' regions are very similar, 
both in shape and intensity (see Fig. ). Only the spectrum 
of the 16" region, although in good agreement with the 
two previous ones above 2 keV, is clearly weaker in the 
soft band. 

3.2.1. EPIC spectrum of the inner region 

The EPIC spectrum of the inner circular region of 
NGC 1808 is displayed in Fig. The extraction region 
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Fig. 4. The Chandra ACIS hard (E>2 keV) and soft (E<1.5 keV) X-ray images of the two nuclear point-like sources. 
The eastern source, SI, is softer than the western one, S2. The nuclear location from 2MASS, indicated by an arrow 
in the images, is closer to source SI. The circles indicate the regions used for spectral extraction (Sect. EO^ . 

Lockman, 1990), a power law with index T = 0.79±£;^ 
and a mekal thermal component with a temperature of 
kT=0.580±g;{£§ keV absorbed by a Hydrogen column of 

Nh=2.01o!4 x 10 21 cm~ 2 . The value of x 2 is unacceptable, 
but an intrinsic absorption to the power law component 
does not improve the fit; the inferred upper limit for the 
Hydrogen column density is 8 x 10 20 cm -2 . The fit only 
improved, in terms of \ 2 , when some abundances were 
allowed to differ from the solar values. The abundances of 
Ne, Mg, Si, and Fe, left free in the fits, are 1.6±J;g, 2.2±^;|, 
2.1+Q'g, and 0.631°'^ times the solar value (Zq), respec- 
tively. Fits with other thermal models such as Raymond- 
Smith or bremsstrahlung emission lead to higher values of 
X 2 - No Fe Ka line was significantly detected. An upper 
limit of 80 eV was calculated for the equivalent width of 
a narrow iron line. Table shows the goodness of the fit 
together with the values and uncertainties for each free 
parameter. 
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Fig. 5. EPIC spectra of the three circular regions and the 
annulus (see text for further details). 



has a radius of 16" (850 pc), centred on the maximum 
peak of the X-ray emission. MOS1 and MOS2 spectra 
and response matrices were combined by merging these 
to maximize the signal-to-noise level. Subsequently, the 
EPIC-pn and combined MOS spectra were fitted simulta- 
neously. Figure shows the 0.35-9 keV spectrum, the best 
fit model with a statistical quality of xt — 0-99 for 174 
degrees of freedom (dof), and the residuals. 

Simple models as a power law, a Raymond-Smith 
(Raymond and Smith, 1976 ), mekal (Mewe, Lemen, & 
van den Oord 1986), and bremsstrahlung components 
provide unacceptable fits. The best fit model includes, 
apart from the fixed Galactic foreground absorption with 
a column density of N H = 3.23 x 10 20 cur 2 (Dickey & 



3.2.2. The Chandra spectra of sources SI and S2 

The ACIS image in Fig. [3] shows a double-peak nuclear 
structure surrounded by a few point-like sources inside 
the inner, r=16", region defined for EPIC and analysed 
in the previous sub-section. We analysed the standard ex- 
tracted spectra obtained with psextract of the two circu- 
lar, r=l'/2, regions SI and S2 (see Fig.0J|. The background 
was extracted from a clean circular region with a radius of 
10" located North-East from the sources in the same 
CCD. For the two regions, models with a single component 
were rejected due to the high \ 2 value obtained. 

The best fit model obtained for the observed spectrum 
of the hard source, S2, includes an absorbed mekal com- 
ponent and an absorbed power law. The values of the 
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Table 1. Parameters and goodness of the best fit models for the r=16", r=62'/5 circular regions, and the 16'.'5-35'' 
annular region. 
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keV 
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keV 
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Ne 1.6+^ 
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Si 2.1±J;g 
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2.0 (f) 


0.79(f) 
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0.58 (f) 


1.1 ± 0.4 


0.096^° 




Ne 0.61S-; 
Mg L0lg:5 
Si 1.2lS'« 


X 2 = 135 for 88 dof 
















Fe 0.37±» : i 2 


6275 (3.3 kpc) 

X 2 = 222 for 222 dof 
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0.79(f) 


6.5 ±0.2 




4.8tM 


H+ u u:i 

u - li -0.02 


4± s » 


Ne 1.5tS'S 

Mg i.9±S:S 

Si 1.9l° ;S 
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The best model for 16" region is phabsGai(zpow±zphabs*vmekal) and the best fit model for 16" 5-35" and 62" 5 regions 
is phabGai(zpow±zphabs*(vmekali±vmekah/). Normalisation in units of: (') photons/keV/cm 2 /s at 1 keV and (*) in 
i.DHil^ Jn e n H dV. 




channel energy (keV) 



Fig.( 

the r= 



EPIC spectrum, 
L6" inner region. 



best fit model, and residuals of 



component with photon index and N# fixed to the val- 
ues obtained for the best fit model of source S2. In or- 
der to model the soft range, we considered an absorbed 
mckal thermal model with free temperature. The result- 
ing fit is unacceptable, xt — 1-9. The addition of an extra 
thermal component with kT=0.087 keV equal to the one 
detected in source S2 does not significantly improve the 
fit. The best fit model is obtained with an unabsorbed 
power law with r = 0.99±g;it and a moderately absorbed, 
N/f = 5.21q'" x 10 21 cm~ 2 , mekal component with a tem- 
perature of kT=0.58ig;|£ keV and the abundances of Ne, 
Mg, Si, and Fe fixed to the values obtained in the fit of the 
EPIC spectrum of the inner region. If an absorption com- 
ponent above the Galactic value is applied to the power 
law, the value of N# for that additional component is com- 
patible with zero. The observed data, the best fit model, 
and the residuals are shown in Fig. 01. 



parameters and the goodness of the fit are presented in 
Table [21 and the observed data, the best fit model, and 
the residuals are shown in Fig. Et>. It is worth noting the 
high values of the absorbing column densities for both 
components, Nh = (2 — 3) x 10 22 cm" 2 . The index of 
the power law found is T = 1.2 ± 0.3, and the temper- 
ature of the thermal component is kT=0.087±Q;o°g keV. 
The thermal component in S2 is significantly colder than 
the one detected in the EPIC spectrum of the inner re- 
gion, with kT=0.58±0.02 keV. The fit does not improve 
by adding an extra thermal component that could account 
for the thermal emission observed in the inner region of 
XMM-Newton, kT=0.58 keV. Neither does it improve by 
varying the abundances from the solar value. 

Figure. 0] shows that above 2 keV, source SI is dim- 
mer than source S2, while below 1.5 keV it is the oppo- 
site. Assuming that the low hard emission of SI could 
be due to contamination of the harder source, S2, we 
tested a model which includes an absorbed power law 



3.2.3. Simultaneous analysis of EPIC-pn and RGS 
spectra 

In order to study the nuclear emission of NGC 1808 in 
greater detail, we took advantage of the high spectral res- 
olution of the RGS data (0.35-2.5 keV). 

Figure [5] compares the combined RGS1 and RGS2 
spectrum of NGC 1808 with that of M82, a proto-typical 
starburst galaxy. Both spectra were generated with the 
SAS task rgsfluxer applying the standard extraction 
techniques. The RGS spectrum of NGC 1808 shows no 
continuum emission above the noise level, but it does ex- 
hibit a number of emission lines. Figurc[S]shows that these 
have wavelengths and relative intensity ratios very similar 
to the strongest lines identified in M 82 (Read and Stevens, 
2002). 

We determined the wavelength of the emission lines ob- 
served with the RGS by fitting Gaussian profiles to them. 
The rcdshift-corrected locations of the fitted Gaussians are 
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Table 2. Parameters and goodness of the best fit models for sources SI and S2 observed with Chandra. 
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Source S2 

x 2 = 41 for 28 dof 
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The best models for regions SI and S2 are phabsGai(zpow+zphabs mc k a i*vmekal) and 
phabsGai(zphabs pw i w *zpow+zphabs me kai*mekal), respectively. Normalisation in units of: (') photons/keV/cm 2 /s at 1 
keV and (*) in 4 ^° (1 + z) -i J n e n H dV. 



SI Chandra Spectrum 



S2 Chandra Spectrum 





channel energy (keV) 



channel energy (keV) 



Fig. 7. Chandra spectra of source SI (a) and source S2 (b), the corresponding best fit models, and the residuals. Note 
that the y-axis scale is the same on both plots to allow comparison of the spectra of both sources. 



less than 0.02 A from the theoretical value of the lines that 
we have identified. The fitted line widths are within the 
range of 0.02-0.04 A, compatible with the instrumental 
resolution of the RGS. The detected lines are thus identi- 
fied as: the Lya emission line from O VIII , transitions of 
He-like NelX ions, and Fe L emission lines from FeXVII 
and Fe XVIII . The weaker lines visible in the M 82 spec- 
trum were not detected in NGC 1808, most probably due 
to the lower S/N. The similarity in the line ratios suggests 
that there is a common origin of the soft X-ray emission in 
both galaxies, i.e. thermal emission from a hot, extended 
gas component, as shown for M82 (Read and Stevens, 
2002). This result confirms the detection of the starburst 
in NGC 1808 in our data, which dominates the soft X-ray 
emission. 

The RGS spectra of NGC 1808 and M82 in Fig.EJshow 
only a qualitative comparison. In particular, in the reduc- 
tion process of the data, calibration effects affecting the 
line profile were not considered, and therefore the wings 
of the lines are not properly determined. The most impor- 
tant effect of this issue is that the flux of the lines could 
be underestimated. The ratio of fluxes of the detected 
lines can be used to derive physical properties of the gas 
where the lines originate. In order to perform an accurate 
study of the properties of the emission lines, we analysed 



the RGS spectrum, taking the response of the detector 
into account and therefore avoiding the uncertainties ex- 
plained before. Using XSPEC v. 11. 3, we simultaneously 
analysed the RGS and the 16" region EPIC-pn spectra. 
Figure [5] shows that the soft X-ray emission of NGC 1808 
can be explained solely by the contribution of emission 
lines. Consequently, we tried to model the pn-RGS spec- 
tra with an absorbed (Nh = lti 7 x 10 19 cm~ 2 ) power law 
(r = l.l±oj and A p i = 9.4+°,'?,xl0- 5 photons/kcV/cm 2 /s 
at 1 keV) which accounts for the hard band emission and a 
number of lines accounting for the soft band emission. The 
emission lines were modeled with Gaussian profiles with 
widths fixed to their instrumental value, both in pn and 
RGS. The wavelengths were also fixed to the theoretical 
values of the most likely identifications in order to reduce 
the degrees of freedom. The best fit model consists of 13 
emission lines and a power law with a value for \ 2 = 307 
for 209 dof. The energy, intensity, and identification of 
each line are shown in Table [21 and the fit in Fig. 

The line ratios provide a diagnostic of the conditions 
of the gas. In particular, the ratios between lines of iron 
in different degrees of ionization degrees allow for the best 
temperature diagnostics. Table [3] gives all the ratios rela- 
tive to the strongest isolated line, i.e. FeXVII (15.01 A). 
We therefore use the observed FeXVII line and compare 
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Fig. 8. RGS fluxed spectra of M82 and NGC 1808, upper and lower panel, respectively. The dotted line corresponds 
to the errors. The identification of M82 lines is from Reads & Stevens (2002). 



it with FeXVIII and others with higher ionization states 
to infer the plasma temperature. We use the web version 
of the ATOMDB v. 1.3, http://cxc.harvard.edu/atomdb 
and conclude that a collisionally ionized plasma with 
kT~0.4 keV (~ 5 x 10 6 K) has a line spectrum compatible 
with the observed spectrum of NGC 1808. This is not sur- 
prising, given our previous results using EPIC data alone. 

3.2.4. EPIC spectrum of the non-nuclear extended 
region 

In order to investigate the non-nuclear extended emission 
of NGC 1808, the EPIC-pn spectrum of an annular region 
was extracted and analysed. The inner and the outer radii, 
16'.'5 (875 pc) and 35" (1.9 kpc), were chosen to minimize 
the contribution of the so-called inner region, i.e. r=16", 
and to exclude the CCD edges. 

For the spectral fit, we have assumed that the emis- 
sion in this region is contaminated by the emission from 
the inner region. Therefore, the fitting model includes a 
power law with an index fixed to the value obtained in the 
fit of the inner region spectrum, i.e. T = 0.79, and an ab- 
sorbed mekal component with Njj fixed to 2xl0 21 cm~ 2 



and kTi fixed to 0.58 keV. The normalisations of both 
components were left free. An excess with respect to the 
model is observed at low energies. This feature can be well 
fitted with an extra mekal component with a temperature 
kT n ~ 0.1 keV. The abundances of Ne, Mg, Si, and Fe 
have been left free in the fits and their values are com- 
patible within the errors with the abundances obtained 
for the inner region. The values of all parameters and the 
goodness of the fit can be found in Table Q Figure E3 
shows the data, the best fit model, and the residuals. 

Comparison of the inner and annular regions shows 
that the annular spectrum includes a component not de- 
tected in the inner region, the low temperature gas, which 
can therefore be associated with a non-nuclear, extended 
emission component. 

The results of the analysis confirm our expectation 
that the spectrum of the annular region is contaminated 
by some emission from the inner region. This is especially 
true for the power law, which is fainter by a factor of 
~6.5 compared to that in the inner region, while the hot- 
ter thermal gas contribution may be partially due to ei- 
ther extended emission or discrete sources in the annulus, 
as its normalisation is only a factor of 3.4 weaker than 
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Energy 


Wavelength 


Intensity 


I/IpleXVII 


Identification 


keV 


A 


ph cm - s _1 






1.99 


6.23 


a° 2+ oo2 


0.01 


SiXIV Lya 6.182 A 


1.87 


6.63 


0.69+° : " 

— U.lo 


0.21 


SiXIII Hea 6.65(6.69)6.74 A 


1.47 


8.43 


0.53 q g 


0.16 


MgXII Lya 8.42 A 


1.34 


9.25 


1.09+°f 9 


0.33 


MgXII Hea9.17(9.23)9.31 A 


1.13 


10.97 


3.04+g : £ 

— U.oo 


0.91 


FeXVII+FeVIII 


1.02 


12.15 


2.34±g'£ 


0.70 


NeX Lya 12.134 A+FeXVII 


0.917 


13.52 


O qo+0.30 
/ - 88 -0 29 


0.86 


NelX Heal3.45,13.55,13.70 A+FeXVII 


0.876 


14.15 


i on+0.26 
i-O^_0.26 


0.39 


FeXVIII 14.21,14.37,14.23 A 


0.823 


15.06 


o oa+0. 30 
0.01_ 29 


1.00 


FeXVII 15.01,15.26 A 


0.769 


16.12 


i qn+0.26 
J "°' fi — 0.26 


0.40 


FeXVIII 16.00,16.07 A+OVIII Ly/3 16.01 A+FeXIX 


0.724 


17.12 




0.75 


FeXVII 16.78,17.05,17.10 A 


0.651 


19.04 


1 4O+0.34 


0.45 


OVIII Lya 18.97 A+CaXVII+CaXVI 


0.565 


21.94 


17+ 042 
u - x ' -0.17 


0.05 


OVII Hea 21.6,21.8,22.1 A 



phabs(zphabs*zpow+13zgauss) xt=lA7 209 dof 




r=16'/5 and r=35" from the centre of the point spread 
function. 
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Fig. 9. High-resolution RGS1+RGS2 and EPIC-pn spec- 
tra in both energy (top plot) and wavelength (bottom 
plot) units, together with the best-fitting model which 
includes a power law and 13 Gaussian profiles and the 
EPIC-pn and RGS-l residuals. 



in the inner region. Indeed, the encircled energy function 
for EPIC-pn predicts that about 20% of the total energy 
emitted by the unresolved source will be imaged between 
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Fig. 10. EPIC-pn spectrum, best fit model, and residuals 
of the annular region. 



3.2.5. Integrated EPIC spectrum of the total galaxy 

A circular region with r = 62'.'5 (3.3 kpc) was analysed. 
The extracted EPIC spectrum contains the emission of the 
whole galaxy and the region coincides with the extraction 
radius used in the R OSAT analysis (Dahlem et al. 1994; 
Junkcs et al. 1995). 

Analogously to the annular region, the components of 
the best fit model ( xt — 1-00 for 222 dof) are a power 
law and two absorbed thermal components. The values of 
the various parameters and the goodness of the fit can be 
found in Table ^ Figure ^] shows the data, the best fit 
model, and the residuals. In the lower panel of Fig. ^2 we 
also show the best fit model and its various additive com- 
ponents in Ef(E). The values of the parameters suggest 
that the emission has the same origin as for the annular 
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region; i.e. the power law and higher temperature com- 
ponent are associated with the nuclear emission, and the 
soft thermal emission is associated with the non-nuclear 
extended emission. 
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Fig. 11. EPIC spectrum, best fit model, and residuals of 
the r=62'/5 (3.3 kpc) region (top plot) and the best fit 
model, the various additive components, and the EPIC- 
pn spectrum in Ef(E) (bottom plot). 



3.3. X-ray fluxes and luminosities 

Various absorbed X-ray fluxes and unabsorbed luminosi- 
ties, measured in the three regions of the XMM-Newton 
analysis and the two regions of the Chandra analysis, are 
collated in Table 01 

For all cases, the thermal components dominate in 
the soft X-ray band, while the power law dominates in 
the 2-10 keV band. The innermost part of the galaxy 
(r=16"=850 pc) is responsible for emission in the hard 
band, while in the soft energy band it contributes only half 
of the total luminosity of the galaxy (r = 62'.'5=3.3 kpc). 

In order to study the long term variability in 
NGC 1808, we compared the fluxes observed with XMM- 
Newton in April 2002 with the ones observed with Chandra 



eight months later. To allow comparison, the spectra of 
the nuclear, the annular, and the inner regions were ex- 
tracted from the Chandra ACIS observation. These regions 
are equivalent to the ones analysed with XMM-Newton. 
The spectra of the inner regions were fitted using the best 
fit model obtained from the XMM-Newton analysis, i.e. 
a power law plus an absorbed mekal component. The in- 
dex of the power law, hydrogen column density, and the 
temperature and abundances of the mekal component of 
the best fit model (x 2 = 170 for 137 dof) are compatible 
with the values obtained with XMM-Newton. Therefore, 
no spectral variability was detected between both obser- 
vations. The spectrum of the annular region was also fit- 
ted using the best fit model obtained with XMM-Newton. 
The values of the relevant parameters of the best fit model 
(x 2 = 45 for 36 dof) are again compatible with the ones 
derived from the XMM-Newton observation. 

In the inner region, the comparison of the XMM- 
Newton and Chandra values reveals no sign of variability 
in the soft energy range Chandra flux is only 5% lower 
than the XMM-Newton flux, well within the statistical er- 
rors and calibration uncertainties of both instruments. In 
the hard band, a decrease of 16% in flux is observed at a 
2.5cr level (see Table EJ). 

The Chandra analysis of the two discrete sources de- 
tected within the r=16" region reveals that source S2 is 
responsible for ^50% of the hard emission from the inner 
Chandra region, while the combined emission of sources 
SI and S2 only explains ^15% of the observed soft emis- 
sion in this inner region. It is worth noting that the up- 
per limit of the intrinsic, i.e. absorption-corrected soft X- 
ray luminosity of S2 obtained with Chandra, is more than 
three orders of magnitude higher than the luminosity of 
the whole galaxy. However, considering the uncertainties 
of this measurement, it is compatible with lower values. 
Moreover, taking the absorbed fluxes into account, no dis- 
crepancy is found, indicating that the intrinsic absorption 
components are responsible for the high luminosity mea- 
sured. Although the statistical test shows that the model 
is acceptable, alternative scenarios cannot be ruled out. 

In the annular region, the Chandra flux measurement 
in the soft band is significantly lower than those of XMM- 
Newton. However, Chandra is less efficient than XMM- 
Newton to detect weak extended emission and probably 
the discrepancy is due to the extended emission not de- 
termined by Chandra. 

We also measured the fluxes and luminosities with 
Chandra within an r=62'.'5 aperture to compare the XMM- 
Newton data with Chandra. The Chandra spectrum does 
not require the complexity needed to describe the XMM- 
Newton spectrum: only a combination of a power law and 
one mekal component is sufficient. The lowest temperature 
mekal, kT~ 0.1 keV, found with XMM-Newton is not de- 
tected in the Chandra spectrum due to the small effective 
area of the instrument at such low energies. The values 
of the parameters of the two components included in the 
best- fit model (x 2 =220 for 200 dof) are compatible with 
the XMM-Newton results. 
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Table 4. NGC 1808 X-ray absorbed fluxes and unabsorbed luminosities 



Region 

radius 


Fo.l-2.4 keV 

10" 13 ergs -1 cm -2 


F2-10 fcV 

10" 13 ergs -1 cm -2 


Lo.1-2.4 keV 

10 39 ergs" 1 


L2-IO keV 

10 39 ergs" 1 


16" (850 pc) 


5.8±i'i 


11.2 ±0.4 


18 ±4 


16.1 ±0.6 


Power law 


1.95 ±0.17 


11.0 ±0.4 


3.2 ± 0.3(~ 20%) 


15.6 ± 0.06 (~ 98%) 


Thermal 


3.8+i'J 




15 ± 4(~ 80%) 


0.33±° [{° (~ 2%) 


16" -35" (875-1900 pc)(f) 






8.2±i; ( i 


2.5±0.5 


Power law 


0.33tg;? a 


i-7itS : iS 


0.5±0.2 (~ 5%) 


2.4±0.6(~ 100%) 


Thermal I 
Thermal II 


1.1±0.02 
0-09^^5 


u.u/o_ 002 


4.3+i'J(~ 40%) 
3.4±1:| (~ 55%) 


0.18±° 2 02 (~ 5%) 


62'.' 5 (3.3 kpc) (f) 


6.8+S:? 


12.1±0.6 


26+S 




Power law 


2.1±1.4 


11.9±0.6 


3.3+ 9 3 (~ 10%) 


16.9+g-J (~ 98%) 


Thermal I 


4.61S1 




15±5 (~ 60%) 


0.38lg; 9 8 (~ 2%) 


Thermal II 


0.18+ 1 ,;?! 




7+1 (~ 30%) 




SI (65 pc) 
Power law 
Thermal 


0.631^*0 
41+ 010 


1 o +u - 3 
1 o +0 - 2 

0.045l°; 2 02 


3.9±i;g 

0.39i«;*3 (~ 10%) 
3.51^ (~ 90%) 


i-4±8:i (~ 95%) 

0.074l 9 ; 93 o 7 (~ 5%) 


S2(65 pc) 
Power law 


n 97 +u.uy 

f _0.18 

20 +0 12 


5.4+£ 8 
5.4t°- 5 


< 6 x 10 4 

< 400 (~ 0%) 


9.2 +u 7 8 

9.4l 9 9 7 (~ 100%) 


Thermal 


n7+ 07 

u - u '-0.03 


3xl0~ 5 


< 6 x 10 4 (~ 100%) 




16" (850 pc) Chandra 


d - d -0.9 








16"-35" (875-1900 pc) Chandra 


5 +u ' 4 

U -°-0.3 


< 0.23 






62" 5 (3.3 kpc) Chandra 


6 l +u ' 4 

0.1_ 1 2 


».X_ 1 2 


26lg 


22.2^ 6 


6275 (3.3 kpc) ROSAT 






~ 140 





(f) The label Thermal / corresponds to the highest temperature mekal component and the Thermal II to the lowest. 



As expected, taking the results for the inner region 
into account, a mild variability in the hard band was ob- 
served between the XMM-Newton and Chandra observa- 
tions, while no significant variation of the flux was mea- 
sured in the soft energy range. This variation is dominated 
by the power law component of the inner region. 



In addition to the newly derived luminosity val- 
ues, Table 0] also lists the luminosity measured with 
ROSAT (Junkes et al. 1995) in this region. The val- 
ues indicate an important decrease in the luminosity 
between the ROSAT, XMM-Newton, and Chandra ob- 
servations. Comparing our results (Lo.s-iofceV ~ 3 x 
10 40 ergs" 1 ) for the r=62"5 region with the ASCA ob- 
servations of February 1994 obtained by Awaki et al. 1996 
(Lo.5-iofceV = 2 x 10 40 ergs -1 ), a small increase in the lu- 
minosity is observed. The comparison with the Ginga re- 
sults (October 1990) indicates a decrease of the luminosity 
from L,2-iokeV = 5 x 10 40 ergs -1 for Ginga to £2-10 keV = 
1.7 x 10 40 ergs -1 for XMM-Newton. Nevertheless, existing 
uncertainties in the relative calibration of various X-ray 
detectors could explain the discrepancy of XMM-Newton 
and Chandra with the ASCA and Ginga values, but the 
luminosity measured by ROSAT is 5 times that observed 
with XMM-Newton and Chandra, too high to be explained 
by calibration uncertainties. 



Short-term variability during the XMM-Newton obser- 
vation was not detected. 



3.4. OM optical/UV data and Ha imagery 

OMUV images of NGC 1808 with the U (350 nm), UVW1 
(291 nm) and UVW2 (212 nm) filters were taken in 
parallel with the X-ray observations. Due to technical 
problems, part of the UVW1 image was lost. All im- 
ages are displayed in Fig. ^] The U and UVW2 im- 
ages show a double-peak structure in the centre of the 
galaxy. The locations of the peaks in UVW2 image are 
a? M (J2000)=05:07:41.97; 5p M (J2000)=-37:30:43.3 and 
a 2 3M (J2000)=05:07:42.62; ^ M (J2000)=-37:30:48.1 with 
uncertainties of ~2". FigurelT2lc shows the central region 
of the UVW2 image with the Chandra contours overlaid. 
The locations of the two point-like sources detected by 
Chandra (at ~4" from the nucleus and aligned in the NW 
and SE direction) coincide with the positions of the UV 
peaks observed in the OM image. However, there is no 
UV emission associated with any of the two bright X-ray 
sources at the very nucleus of the galaxy. 

In Fig. ED we have overlaid the OM 212 nm image of 
NGC 1808 on the EPIC 0.5-2.0 keV frame from Fig.HJand 
on an extremely sensitive Ha image. An RGB composition 
of the images in optical (Ha), UV {OM UVW2) and X- 
rays (EPIC) is presented in Fig. ITU 

The correspondence of UV and Ha emission maxima 
indicates that the OMdata trace the massive stars ionizing 
the gas in HII regions in NGC 1808. The Ha image also 
shows several knots located in an outer spiral arm. SF 
is traced in the circumnuclear starburst but also further 
out in the galaxy disk. Similarly, the soft X-ray emission 
observed in the same regions is expected to come from hot 
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Fig. 12. OMUV images of NGC 1808. From left to right and from top to bottom: U filter (350 nm), UVW1 filter (291 
nm), UVW2 (212 nm) images, and the UVW2 central region of NGC 1808 with the overlay of the Chandra contours. 
The location of the nucleus from 2MASS is marked with a cross. The elliptical ring eastwards of the galaxy in U and 
UVW1 images is an artefact due to the light scattered within the OM detector. All images are displayed in logarithmic 
scale with a threshold of 3<r of the mean background value. 



gas in or near these SF regions. Optical/UV emission also 
tracks the bar. 



sion of an AGN (e.g. Junkes et al. 1995), thus suggesting 
the predominance of the starburst. 



4. Discussion 

ASCA and Ginga observations in the past missed the star- 
burst's soft X-ray emission, favouring an interpretation in 
terms of an AGN (Awaki et al. 1996, Awaki & Koyama 
1993), as suggested first by Veron-Cetty and Veron (1985). 
Similarly, the IR emission observed by ISO is dominated 
by this AGN (Sicbcnmorgen et al. 2001). 

ROSAT, on the other hand, detected the soft emission 
from the starburst, being less sensitive to the hard emis- 



Based on the data presented above, the most likely in- 
terpretation is that NGC 1808 hosts both, a circumnuclear 
starburst and an unresolved nuclear X-ray source. Both 
in the imaging (Figs-DEl andEJ) and spectral (Figs. Eld 
and |3J) domains, XMM-Newton and Chandra detect the 
presence of the starburst and the hard unresolved nuclear 
source. In the following we discuss the physical origin of 
the emission of both types of activity: the starburst and 
the unresolved X-ray source. 




Fig. 13. Overlay of our UVW2 (212 nm) filter OM contours on top the EPIC 0.5-2.0 keV frame (left) and an Ha 
image (right). (Contour levels are 2, 7, 9, 10, 15, 20, 30, and 80 counts.) All images are displayed in logarithmic scale 
with a threshold of 3er of the mean background value. 



4.1. X-ray emission of the inner region 
4.1.1. Hard component 

Spectral characteristics 

The slope of the power law fitted to the 16" re- 
gion, r = 0.79^Qgg, is flatter in comparison to results 
obtained for radio-quiet AGN and LINERs. Works by 
Georgantopoulos et al. (2002) and Terashima et al. (2002) 
find values for the power law spectral indices ranging 
from 1.7 to 2.3 for these objects. However, there are also 
objects with lower spectral indices (e.g. PKS 2251+113, 
with T = 0.95 ±0.24, Reeves & Turner, 2001). Analysis of 
the Chandra data reveals that this hard emission is dom- 
inated by one of the two detected nuclear sources, S2 in 
Fig. 21 which is harder. In addition, source S2 is absorbed 
by a column density of N# = 3.1+q' 8 , x 10 22 cm -2 , com- 
patible with values observed in Compton-thin Seyfert 2 
galaxies. This fact suggests that in case one of the nuclear 
peaks observed with Chandra is associated with an 
AGN, the most likely possibility is that it is source S2. 
Interestingly, the best-fit model of the XMM-Newton 16" 
region spectrum does not require absorption in excess of 
the Galactic value. However, this XMM-Newton region 
includes several emitting sources that could blur the 
presence of an absorbing component in the best fit model. 

Luminosity 

The total luminosity of the central region 
with a radius of 16"= 850 pc of NGC 1808 is 
L 2 -iokcV = (1-61 ± 0.06) x 10 40 erg s" 1 . The lumi- 
nosity in the same energy band associated with the 



harder peak with a radius of 1"2=65 pc observed with 
Chandra, i.e. source S2, is L 2 -iokcV = 0.92to;? 8 X 10 40 erg 
s , which is ~ 60% of the emission of the XMM-Newton 
inner region. These values are in the range of Low 
Luminosity AGN (LLAGN) and LINERS (Terashima 
et al. 2002, Georgantopoulos et al. 2002). Moreover, in 
both cases the luminosity associated with the power 
law component is in good agreement with the typical 
values of low luminosity Seyfert galaxies and LINERs. 
Terashima et al. (2002) found in a sample of 53 objects 
of these types observed by ASCA that the luminosity 
varies within the range L2-wkeV = 5 x 10 39-40 ergs -1 . 
If only LINERS are considered, the mean luminosity 
obtained is L 2 -wkeV = 7 x 10 39 ergs -1 (Ho et al. 2001), 
which is lower than the one derived for NGC 1808. In the 
soft band, a sample of LINERS observed with ROSAT 
(Komossa et al. 1999) shows luminosities in the range 
of Lo.24-2fceV = 10 38-41 ergs -1 , in agreement with 
luminosity measured for the NGC 1808. 

Variability 

We studied the variability of NGC 1808. While no 
short-term variation of flux was detected during the 
XMM-Newton observation, medium- and long-term vari- 
ability has been suggested in both soft and hard bands. 

NGC 1808 was observed by XMM-Newton and 
Chandra with an interval of eight months. The measure- 
ments indicate the soft X-ray flux remained constant, tak- 
ing the uncertainties into account. A decrease of 16% (at 
a 2.5(7 level) in the measured flux is observed in the hard 
X-ray band. In particular, the power law component of the 
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innermost region spectrum of NGC 1808 is responsible for 
this variability. 

The long- and medium-term flux variations measured 
over the years by various satellites are most probably 
caused by variability of the nuclear source. In the soft 
energy band, a significant decrease of a factor of five has 
been observed when the ROSAT and X MM- Newton (or 
Chandra) luminosities are compared. It should be noted 
that the ROSAT spectral fit is particularly uncertain be- 



cause of the high absorbing column density, which leaves 
a very narrow effective bandpass for the X-ray emission, 
from about 0.5 to 2.4 keV. 

In summary, mild medium-term variability of the hard 
band luminosity was measured between the XMM-Newton 
and Chandra observations. This variation can be at- 
tributed to an AGN or a ULX. Furthermore, the uncer- 
tainties introduced by comparing different satellites does 
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not allow us to firmly conclude that other changes in flux 
are intrinsically due to the nucleus of NGC 1808. 
Ultraluminous X-ray sources 

The study of the hard X-rays emanating from the in- 
ner region indicates that the nuclear X-ray emission itself 
could be due to a LLAGN or a LINER. The observed 
hard band X-ray luminosity, L2-10 fceV = (L61 ± 0.06) x 
10 40 erg s _1 , is very difficult to explain by standard stel- 
lar processes. Nevertheless, non-nuclear point-like sources 
with X-ray luminosities on the order of 10 38-41 erg s _1 
have been detected in nearby galaxies. An alternative ori- 
gin to an AGN nucleus of NGC 1808 could be one of the 
so-called ultraluminous X-ray sources (ULX). 

The two central sources SI and S2 were studied in de- 
tail in Sect. 13.2.21 The nuclear location from 2MASS is 
compatible with either of them being the nucleus. The lu- 
minosity of source S2 of L 2 _io keV = 9.2+°' 8 x 10 39 erg s _1 
is compatible with either a LLAGN or a ULX. Strickland 
et al. (2001) and Roberts et al. (2002) show that the anal- 
ysis of Chandra ULX spectra favours simple power laws 
with indices in the range of 1.8-2.9, higher than the value 
obtained in our analysis of 0.79^g'Qg. Interestingly, one of 
the compact radio sources reported by Saika et al. (1990) 
coincides with the X-ray point-like source SI. 

Although more data analysis is necessary to reach any 
firm conclusion, this finding leaves open the possibility 
that the nucleus of NGC 1808 could be a ULX. 

4.1.2. Soft component 

Emission distribution 

The soft X-ray emission of the inner region of 
NGC 1808 is extended (see e.g. Fig. HJ. The elongated 
emission is detected at energies up to 4.5 keV. Moreover, 
the Chandra images clearly show the simultaneous pres- 
ence of a double emission peak and diffuse, or at least 
unresolved, X-ray emission around the two peaks. The ex- 
tended emission is likely to have been created by star for- 
mation processes. Our XMM-Newton data permit a clear 
spectral distinction between the emission associated with 
the starburst and the additional nuclear sources in the 
r=16" (850 pc) region, see Fig. The starburst, associ- 
ated with the mekal component in the spectral fit, dom- 
inates the soft energy band (below 1 keV) , see Tables ^ 
and 01 The starburst also shows up in the form of radio 
continuum and IR knots (Saikia et al. 1990, Kotilainen et 
al. 1996), bright Ha emission, and bright circumnuclear 
UV continuum from massive stars, as displayed in Fig. 1131 
In particular, the position of the 2MASS nucleus coincides 
perfectly with one of the radio compact sources detected 
by Saikia et al (1990), explained as SNR or complexes of 
unresolved SNR. The nucleus is prominent in X-rays and 
in the optical regime (see Figs. HI and 1 121 a). 

However, while the OM UVW2 image shows that the 
UV emission close to the position of the nucleus is very 
weak (see Fig. El c )- These facts suggest that the nucleus 
is highly absorbed. The spectral analysis of source S2 



shows a high value of the equivalent hydrogen column, 
N# ~ 10 22 cm~ 2 , which can explain this absence of UV 
and soft X-ray emission associated with the hard X-ray 
peak. 

Plasma temperature 

The best fit to the EPIC data contains a ther- 
mal plasma component with a temperature of kT 
= 0.580± ;^§ keV. This corresponds roughly to the 
"medium" component identified in the spectra of other 
starburst galaxies (Dahlem et al. 1998; see also Weaver 
et al. 2000 and Dahlem et al. 2000). It is consistent 
with 0.5 ± 0.2 keV (Junkes et al. 1995), which was 
measured from ROSAT data, based only on about 600 
photons. It is interesting to note that the surprisingly 
high Hydrogen column density, on the order of 8 x 10 21 
cm~ 2 , in the best-fitting model to the ROSAT data is 
also corroborated by our results. 

X-ray luminosity and SFR 

The luminosity of the starburst component, which we 
identify with the thermal spectral components, is domi- 
nated by the emission in the soft energy band. The ob- 
served X-ray luminosity, £0.2-4 fceV ~ 1.3 X 10 40 erg s _1 , 
which represents ~70% of the emission associated with 
the starburst component of the whole galaxy, leads to a 
far-infrared-to-X-ray luminosity ratio of log(Lx/^FiR.) = 
-3.8 2 . This value is consistent with those obtained by 
Heckman et al. (1990) for a sample of six starburst galax- 
ies, which are in the range between -3.7 and -4.3, and also 
with the results in Mas-Hesse et al. (1995), who obtained 
a mean value of -3.33 for a sample of starforming galaxies. 

Ranalli, Comastri & Setti (2003) deduced that the X- 
ray luminosity is a tracer of SFR, based on a sample 
of nearby galaxies. As was pointed out by Persic et al. 
(2004b), the luminosity in the 2-10 keV band is not a pre- 
cise SFR indicator, due to the contamination by low-mass 
X-ray binaries (LMXB). These objects are bright in the 
X-ray band over long timescales, tx ~ 10 7 yr, and there- 
fore do not trace the instantaneous SFR. For this reason, 
we use the soft X-ray luminosity as a tracer of the SRF, 
with the SFR-L0.5-2 keV relation in Ranalli, Comastri & 
Setti (2003): 

SFRxiMayr- 1 ) = 2.2 x 10- 40 L .5-2 keV (1) 

The Lo.5_2 keV associated with the starburst in the 
r=16"(850 pc) region, as derived from the mekal compo- 
nent contribution, is on the order of LQf_ 2fceV - ~ 1 x 10 40 
erg s _1 and the SFRx, therefore, is ~ 2.5 M Q yr _1 . 

Metallicities 

We detect several key emission lines, which for the first 
time provide us with a good measure of the abundances of 
some elements fSect. r?~2.1|) . The RGS spectrum shows a 

2 Lfir = 2tt D 2 FIR = 9.3 x 10 43 erg s" 1 (Dahlem et al. 
1992), where FIR = 1.26 (2.58 f 60 + /100), ho being the IRAS 
60 /im flux density and /100 the 100 /im flux density. 
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very weak continuum emission with several emission lines. 
The values obtained indicate slightly super-solar metallic- 
ities for Ne, Mg, and Si and solar or sub-solar abundance 
for Fe, which is consistent with the results by Dahlem et 
al. (1998) and Weaver at al (2000). 

In this context the absence of a thermal Fe Ka line 
at ^6.5 keV energy is noteworthy. From our data in the 
annular region, we derive a 3-cr upper limit on the equiv- 
alent width of 170 eV and on the flux of 4 x 10~ 7 pho- 
tons cm, _2 s _1 . This result is compatible with any possible 
iron line emission entirely due to the mekal component. 
We have used ATOMDB v. 1.3 to derive the expected flux 
of the emission lines for a gas with temperature and nor- 
malisation as obtained from the single mekal component 
fitted to the annular spectrum, i.e. kT=0.53lg Qg keV and 
normalisation A mc ]< a i = 2.2 x 10~ 4 . The total flux of the 
most prominent iron lines (FeXXI-6.505 keV and FeXXII- 
6.504 & 6.57 keV ) is - 2 x 10" 4 photons cm^s' 1 (with 
an emissivity of ~ 5 x 10~ 22 photons cm 3 s _1 ), consistent 
with the upper limit measured. 

4.2. X-ray emission outside the central region 

In Sect. 13.21 we show the results of the spectral analysis of 
the 875 pc-1.9 kpc annular region. The best fit model in- 
dicates that, apart from the contribution of the nuclear 
emission, the annular region includes an extra thermal 
component which can be explained by a mekal model with 
kT=0.10+Jj;[| 2 keV. The analysis of the total galaxy X- 
ray emission also indicates the presence of a soft thermal 
model with a similar temperature, kT=0.11^Q Q2 keV. 

Although weak, the contribution in the soft band of 
this component (Lo.i-2.4feeV = 7.7 x 10 39 erg s _1 ) to the 
soft X-ray luminosity of the whole galaxy is significant, 
representing around 30% of the total soft X-ray emission 
of NGC 1808. In the hard band, the contribution of this 
low temperature thermal component is negligible. Using 
Eq. 1 (see Sect. 14.1. "2*1 for further details), the observed X- 
ray luminosity associated with both mekal components of 
the annular region, Lo.5_2 keV ~ 4x 10 39 , can be translated 
into a SFR of -1 M Q yr" 1 . 

Far-infrared emission is also a good tracer of the 
star formation; using the relationship by Kennicutt et al. 
(1998), we derive for NGC 1808: 

SFR FIR (M Q yr- 1 ) = -4.5 x 10- 44 L FIR (2) 
e 

where e is the fraction of UV/optical flux emitted by 
the stars and absorbed by the dust and remitted in 
the IR, which is close to unity even for low values of 
reddening (Mas-Hesse & Kunth 1991). Assuming that 
the whole FIR emission is due to the starburst, we ob- 
tain SFRp-jjj ~4 M^Qyr" 1 . This result is consistent with 
the one obtained through the X-ray emission, combin- 
ing the values derived from the nuclear region (SFRx ~ 
2.5 Mq yr -1 , see Sect. 14. and from the annular region 
(SFRx ~ 1 Af Q yr _1 ). It suggests that the contribution 
of the AGN to the FIR emission is small. 



Table 5. Parameters for the hot gas component in the 
extra-nuclear starburst region of NGC 1808 



(cm' 3 ) (M Q ) (yr) 

3.0 x W~ 2 /y/rj 6.3 x 10 6 ^/7j 6.3 x lO 7 ^ 



The second, soft mekal component detected in both the 
annular and the r=3.3 kpc region is most likely thermal 
emission associated with the halo. Such a thermal compo- 
nent (kT=0.2-0.3 keV) is found in all starburst galaxies 
in the sample by Dahlem et al. (1998) and Weaver et al. 
(2000) and also by us in other starburst galaxies observed 
with XMM-Newton (e.g. Dahlem et al. 2003 and Ehle & 
Dahlem 2004). The presence of both, the prominent ver- 
tical dust filaments in NGC 1808 and the gaseous outflow 
inferred by the kinematics of HI gas (Koribalski et al. 
1993), suggested that this component would most likely 
be present in NGC 1808. 

If the diffuse X-ray emission in the soft energy band 
is assumed to be due to thermal cooling of a hot gas in 
collisional ionization equilibrium 3 , it is possible to esti- 
mate the gas density n e , mass m gas , and cooling time t of 
that plasma. To this end we make use of a formula given 
by Nulsen et al. (1984), i.e. L x (soft) = 0.812A(T)n^Vr]. 
The unknown filling factor rj allows for some dumpiness 
of the gas filling the emitting volume V, and was found 
in simulations (Avillez & Breitschwerdt, 2004 and refer- 
ences therein) to be in the range of 0.2-0.4, only modestly 
increasing with higher SN rate. 

As explained in Sect. 3.2.2, the spectrum of the non- 
nuclear extended X-ray emission is complex and a best fit 
was obtained with a model consisting of a foreground ab- 
sorbed power law plus (at least) two internally absorbed 
mekal components, corroborating the fact that the hot 
gas is a multi-temperature mixture. Nevertheless for the 
purpose of estimating hot gas parameters, we fitted this 
emission component with a single 'characteristic' temper- 
ature that was found to be 0.53±g;g| keV. For such a gas 
temperature of the starburst component (~ 6.2 x 10 6 K) 
Raymond et al. (1976) give a cooling coefficient A(T) of 
~ 5 x 10 -23 erg cm 3 s _1 . For the unabsorbed luminosity of 
the hot gas component we make use of the total ( Thermal 
I + IT) Lo.i-2.4 keV = 7.7 x 10 39 erg s" 1 (see Table EJ). 
Assuming for the emitting volume a torus geometry with 
an inner radius of 875 pc and outer radius of 1.9 kpc, the 
calculated gas parameters are presented in Table El 

The mass of the detected hot gas is only ^ 0.5% of 
the total gas mass 4 . The cooling time, r, is comparable 
to typical duty cycles of starbursts, which are estimated 

3 See, however, the discussion in the validity of this assupn- 
tion in Breitschwerdt & Schmutzler (1999) 

4 The total gas mass calculated using H/ observations is > 
1.8 x 1O 9 M (Dahlem et al. 2001). 
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to be on the order of a few times ~ 10 7 yr (Rieke ct al. 
1988). 

No significant X-ray emission was found, at the sensi- 
tivity of our current data, in the intergalactic space around 
NGC 1808. Similarly, no HI and Ha emission was found 
by us (Dahlem et al. 2001). 

5. Summary and conclusions 

The work presented in this paper shows evidence of the 
presence of starburst activity and a hard unresolved source 
within the inner r=850 pc region in NGC 1808. Analysis 
of the optical-to-X-ray emission proves the co-existence 
of thermal diffuse plasma, non nuclear point-like sources 
and, of a LLAGN or a ULX. 

The XMM-Newton EPIC images show extended X-ray 
emission for energies below 4.5 keV. The elongated shape 
observed in the X-ray band follows the same orientation as 
the Ha and optical-UV emission, suggesting that SF con- 
tributes significantly to the soft X-ray emission. At least 
four non-nuclear point-like sources were detected in the 
Chandra ACIS image of the inner 250 pc of NGC 1808. 
The high spatial resolution Chandra image allows the de- 
tection of a double-peak nuclear structure. The locations 
of both sources are compatible with the position of the 
nucleus as derived from 2MASS. The spectral character- 
istics of these sources show that one is dominated by hot 
gas emission, indicating starburst activity, while the other 
is a highly absorbed hard X-ray source. 

The spectral analysis of the XMM-Newton data com- 
pletes the scenario outlined by X-ray imaging. Comparing 
EPIC-pn spectra of several regions in the central part of 
NGC 1808, it is inferred that the bulk of the emission orig- 
inates from the nucleus, although in the soft X-ray band 
the contribution of the circumnuclear regions is not neg- 
ligible. 

The X-ray EPIC spectrum of the inner nuclear re- 
gion (r=16", 850 pc) is explained by a power law, T = 
0.791q'qq, accounting for the hard X-ray emission and as- 
sociated with an unresolved point-like source. The ab- 
sorption measured in the Chandra spectral analysis of 
source S2 is compatible with typical values of Compton- 
thin Scyfert 2 galaxies. However, this absorbing material 
docs not show up in the XMM-Newton spectrum, prob- 
ably because it includes contributions from different re- 
gions, making it difficult to separate absorption and soft 
emission components. 

In order to explain the soft X-ray emission of the in- 
ner region, it is necessary to include emission due to a 
diffuse plasma with kT=0.580±£;^ keV. The data ob- 
tained with RGS exhibit emission lines similar in wave- 
length and relative intensity ratios to the ones found for 
the prototypical starburst galaxy M82. The temperature, 
the luminosity, and the various ratios of emission lines are 
consistent with being originated in a collisionally ionized 
plasma associated with SF regions, which dominates the 
total emission spectrum including the continuum in the 
soft X-ray regime. 



The luminosities measured with XMM-Newton for the 
r=850 pc inner region (L 2 -iofceV = (1.61 ± 0.06) x 
10 40 erg s _1 ) and with Chandra for source S2 (L2-10 keV = 
0.92±° ° 7 8 x 10 40 erg s" 1 ) are both in the range of LLAGNs 
but are also compatible with the values of ULX sources. 
A decrease in the 2-10 keV flux on the order of 15% in an 
interval of eight months was measured when comparing 
XMM-Newton and Chandra observations, while no change 
in flux was detected in the soft energy band. Both the 
luminosities and the variability detected are compatible 
with the unresolved source being a LLAGN or a ULX. 

Halo emission of NGC 1808 could also be de- 
tected: a softer thermal component, kT— 0.1 keV, with 
L .i-2.4fcev ~7x 10 39 erg s" 1 (-30% of the total soft 
X-ray emission) is necessary to explain the spectrum at 
distances larger than 875 pc from the centre. 

The multiwavelength analysis of the nuclear region 
of NGC 1808 performed in this work allows us to detect 
and for the first time to unambiguously disentangle the 
contributions of an unresolved nuclear X-ray source and 
the starburst regions within the r=850 pc region. Further 
analysis of better quality data is required to determine the 
exact nature of the unresolved nuclear source. 
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